Propagation of transverse intensity correlations of a two-photon state 
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The propagation of transverse spatial correlations of photon pairs through arbitrary first-order 
linear optical systems is studied experimentally and theoretically using the fractional Fourier trans- 
form. Highly-correlated photon pairs in an EPR-like state are produced by spontaneous parametric 
down-conversion and subject to optical fractional Fourier transform systems. It is shown that the 
joint detection probability can display either correlation, anti-correlation, or no correlation, depend- 
ing on the sum of the orders a and /3 of the transforms of the down-converted photons. We present 
analytical results for the propagation of the perfectly correlated EPR state, and numerical results 
for the propagation of the two-photon state produced from parametric down-conversion. We find 
good agreement between theory and experiment. 



I. INTRODUCTION 

The discussion about non-local correlations between 
properties of two separated particles began in part with 
the famous EPR paper in which Einstein, Podol- 
sky and Rosen showed that the position and momen- 
tum of two correlated particles could be used to con- 
struct a paradox between quantum theory and intuitive 
concepts like locality and the reality of physical proper- 
ties. Continuous variable (CV) entangled states similar 
to the EPR state appear in a number of physical systems, 
including field-quadrature correlations of two modes of 
the electromagnetic field [l, 0, 0, Hi j spatial variables of 
pairs of photons [1, 0|, and others [1, i, EH- This 
has allowed for experimental realization of the original 
gedanken experiment proposed by EPR 0, @, CV 
entanglement of the EPR type has been shown to be 
useful for a number of quantum information tasks p^ . 
One benefit to the study and use of CV's is access to a 
Hilbert space of larger dimension, which is advantageous 
for quantum cryptography p^ . IT4 | and fundamental tests 
of quantum mechanics [l5|. 

EPR-like spatial correlations can be identified by the 
violation of the inequality [l^, 



A2(pi|p2)A2(qi|g2) > 



4' 



(1) 
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where A'^{pi\p2) represents the uncertainty in variable 
pi of system 1 conditioned upon measurement of sys- 
tem 2 at p2- A^(pi|p2) is the variance of the conditional 
probability P{pi\p2) for a fixed value of p2- Similarly, 
^'^{Qi\Q2) is the variance of the conditional probability 
P{qi\q2), where qi and (72 are the Fourier conjugate vari- 
ables of pi and p2- If inequality ^ is violated, one could 
infer cither pi or qi from conditional measurements p2 
or q2, with less uncertainty than the Heiscnberg uncer- 
tainty principle would allow. In recent experiments [6|, 0] , 
measurements of this type were performed in the coinci- 
dence counting regime using photons from spontaneous 
parametric down-conversion (SPDC). The transverse po- 
sition and momentum were determined by measuring the 
intensity distributions in the near and far field, respec- 
tively. As inequality ([T]) deals with EPR non-locality, it is 
generally more restrictive than those involving variances 
of center of mass and relative variables, which identify 
non-separability of continuous variable systems (isl . IToj . 

The spatial correlations of photon pairs produced by 
SPDC present a rich playground to investigate CV cor- 
relations with relatively simple linear optical systems 
[23,[Ill,|23,|2l,[2l,[2i,[26[. In SPDC sources, photon pairs 
generally display an intensity correlation in the near field 
(source), due to the localized emission of the photon pair: 
the photons are "born" from the same pump photon, so 
both photons are detected at nearly the same position 
in the source plane. As the entangled two-photon state 
propagates, this spatial correlation evolves to an anti- 
correlation in the far field. Consequently, if photon 1 is 
detected at position p in the far-field, photon 2 will be 
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found near —p. The far- field anti-correlations are due 
to the phase matching (momentum conservation) in the 
non-linear SPDC interaction. The spatial correlation in 
the near-field and anti-correlation in the far-field have 
been previously observed in Ref. [23| ■ The switch from a 
near-field correlation to a far-field anti-correlation raises 
the question as to what type of correlation is present 
at intermediate distances in between the near and far- 
field regions. Recently, Chan et al. [1^ showed that the 
correlations can "migrate" entirely to the phase of the 
two-photon wave function, and consequently the condi- 
tional intensity distribution may display no correlation 
at all. In Ref. [1^, it was shown that it is always pos- 
sible to detect transverse entanglement performing only 
intensity correlation measurements, when an arbitrary 
propagation is applied to each of the entangled photons. 

The propagation of the transverse spatial structure of 
an optical field can be accurately described by the Frac- 
tional Fourier Transform (FRFT) This is true for 
any first-order linear optical system. That includes free- 
space propagation alone [sil, [32, [s^l , and also optical sys- 
tems consisting of lenses and free space [SO,, J4] , provided 
that one chooses the appropriate scaling of the transverse 
coordinates. The FRFT is parameterized by an angle a, 
so that a = corresponds to an identity operation and 
a = 7r/2 is the usual Fourier transform. With proper 
scaling of the coordinates, the FRFT is additive, so that 
consecutive FRFT's !Fa and !Ff} can be written as Ta+p- 
This allows one to associate an overall FRFT to an arbi- 
trary first-order linear optical system. 

In the present work, we study the transverse EPR 
correlations of propagating SPDC photon pairs using 
the FRFT. We show theoretically and experimentally 
that the presence of EPR intensity correlation, anti- 
correlation or no correlation depends on the sum of the 
orders a and {3 of the applied FRFT transforms in each 
of the down converted photons. In this way, it is possi- 
ble to engineer the spatial intensity correlations through 
the application of optical FRFTs to the entangled down- 
converted photons. The FRFT describes a canonical ro- 
tation in phase space, and applies to any pair of conjugate 
variables, such as time- frequency [sot or field quadratures 
p^ . Thus, the conclusions drawn here are also relevant 
to other physical systems. 

In section |lll we review the connection between the 
Hilbert space associated to the spatial variables of a sin- 
gle and two-photon field and the Hilbert spaces of point 
particles with two degrees of freedom. This allows us 
to apply the usual quantum formalism for point parti- 
cles in the description of the spatial properties of single 
and two-photon states. In section III A[ we discuss the 
propagation of photons through first-order linear optical 
systems and the use of the FRFT in this description. Sec- 
tion mil introduces the type of two photon state typical 
of the SPDC process, and discusses the propagation of 
transverse correlations under FRFT operations. In sec- 
tion [IV] we present an experiment and results which are 
well described by the theoretical results presented in sec- 



tion [TTTl Finally, we provide some concluding remarks in 
section |Vl 

II. SINGLE AND TWO-PHOTON STATES 

Here wc focus on the spatial structure of a single or 
two-photon field. Thus, for simplicity, wc will assume 
that the fields are paraxial, monochromatic and have well 
defined polarization. The Hilbert space Tii describing the 
transverse spatial degrees of freedom of a single-photon 
state is spanned by the basis {\p) = a){p)\Q)}, where 
|0) is the vacuum state. An arbitrary pure state is then 

\i^) = J dpw{p)\p), (2) 

where p = (px,Py) is the transverse position and w(p) 
is the transverse wavefunction or detection amplitude. 
The basis states {\p)} correspond in second quantization 
to unnormalized states of one photon at position p. It is 
possible to establish an isomorphism between Hi and the 
Hilbert space spanned by position eigenstates of a two- 
dimensional position operator p = (p^, Py) if one specifies 
the action of this operator on the basis states as: p\p) = 

m- 

Alternatively, Tii is spanned by the basis {|q) = 
a^(q)|0)}, where 

aHu)^ I dpe^^-^a^p) (3) 

and q = {qx,(jy) are the transverse components of the 
wave vector k. In this basis the wavefunction u(q) is the 
angular spectrum of the photon field, and is obtained 
by a Fourier transform of the detection amplitude w(p). 
Again, it is posible to establish an isomorphism between 
Hi and the space spanned by momentum eigenstates of 
a two-dimensional momentum operator q= (q^, %) if the 
action of this operator on the basis states is: q|q) = q|q). 
Because the two bases {|p)} and {|q)} are related via a 
Fourier transform similar to the one in Eq. ^ , the posi- 
tion and momentum operators satisfy the canonical com- 
mutation relations [pk,qi] = i6k,i^, where k,l — {x,y}. 
Thus, at the level of quantum kinematics, there is an 
isomorphism between the Hilbert space corresponding 
to transverse spatial degrees of freedom of single-photon 
states and the Hilbert space of quantum states of a point 
particle with two degrees of freedom. The equivalence 
between the classical paraxial wave optics and the non- 
relativistic quantum mechanics of two-dimensional point 
particles is well known [sl] , and also allows one to estab- 
lish the isomorphism at the level of quantum dynamics. 
In fact, for paraxial propagation of the photons along an 
optical axis z, the wave equation that governs the evolu- 
tion of the wavefunction w{p) = {p\w), is a time depen- 
dent Schrodinger equation where the length variable z 
plays the role of time and the wavelength. A, of the pho- 
tons plays the role of Planck's constant [4J| . The analogy 



3 



between paraxial wave propagation and non-rclativistic 
quantum mechanics of a point particle has been well ex- 
plored M, M, [33, m, m . 

The Hilbert space describing the transverse spatial de- 
grees of freedom of two-photon states is simply the ten- 
sor product Hi (E) 7^2 between the Hilbert spaces of one- 
photon states. Thus, Tii®T-l2 is isomorphic to the Hilbert 
space associated to two distinguishable point particles, 
each one with two degrees of freedom. We assume that 
the photons are distinguishable since in principle they 
could be distinguished by their longitudinal direction of 
propagation or their polarization. Therefore, an arbi- 
trary two-photon pure state can be written as 



l^)= jjdpidp2^{pi,P2)\pl)i\p2) 



(4) 



where ^'(pi,p2) = (pi,/02|^) is the normahzed wave 
function and \pi)i and 1/32)2 ^"^^ position eigenstates for 
photons 1 and 2, respectively. Here it is assumed that 
the paraxial approximation has been applied along two 
distinct z axes, one for each single-photon field. 



A. Propagation as a Fractional Fourier Transform 

The most common optical systems are first-order lin- 
ear systems (also called quadratic-phase systems), which 
are composed essentially of sections of free space and thin 
lenses centered on the propagation (2) axis (30l |. Parax- 
ial propagation in these systems is particularly simple: 
the paraxial wave equation corresponds to a Schrodinger 
equation associated with a quadratic hamiltonian, so the 
evolution of the phase space operators are simply given 
by {PAV U^{pA)'^U = M{p,q)'^ (where U is the 
evolution operator associated with the quadratic hamil- 
tonian, and T means transposition). The symplcctic ma- 
trix M is the ray matrix that stems from geometrical 
optics applied to the system. For example, in the case 
of only free propagation the evolution is associated with 
the hamiltonian of a free particle and the matrix M rep- 
resent a linear canonical transformation that correspond 
to a shear in the direction of the transverse momentum 

A great simplification and systematization in the de- 
scription of evolution through first-order optical systems 
is gained by using dimensionless variables p = p/s and 
q = sq, where the real number s has the dimension of 
a length, and is generally a function of the properties of 
the physical system. In this case, free space propagation 
can be described in the paraxial approximation with the 
help of the FRFT [13, S [13. This is due to the fact 
that the paraxial Fresnel diffraction integral, which re- 
lates the light signal between two transverse planes in 
free space, can be expressed using a FRFT if we use 
dimensionless coordinates. The more general case oc- 
curs when we choose different parameters s at the input 
and the output planes. However, in order to identify the 
tranvcrse position and momentum coordinates at these 



planes as belonging to the same phase space, one must 
use the same parameter s. FIG. [l]a) illustrates identifi- 
cation of the FRFT with propagation through free space. 
The diffraction of light from a spherical cap emitter with 
radius of curvature = — i? < to a spherical cap re- 
ceiver with radius of curvature = i? > at a distance 
z from the emitter can be expressed as (32l . Issj 



0,(p) = e-*"/2.F„[0e(p)], 



(5) 



where 4>e{p) = exp{—ik\sp\'^ /2Re)ipe{p) and (f>r{p) = 
exp{—ik\sp\'^ /2Rr)(pr{p) {k = |k| = 2tt/X). (pe{p) and 
fr{p) are the wavcfunctions at planes of observations 
tangent to the emitter's and receiver's spherical caps at 
its vertex point. Here we call the angle < a < tt 
the order of the FRFT. This order a = a{R, z) and the 
adimensionalization parameter s = s{R, z) can be cal- 
culated from the relations s = \J zjk (1 — 5^)"^^'' and 
g= \ — z / R = cos a. Alternatively, given the parameter 
s and the distance z we can estimate the FRFT's order 
a = a(s,z) and the radius of curvature R = R{s,z). 
It is important to note that the quadratic phase factors 
that maps the wavcfunctions at the spherical caps to the 
wavcfunctions at their tangent planes are not important 
if we are concerned only with intensity measurements at 
these planes. 

The description of the propagation of photons through 
first-order optical systems with the help of the FRFT is 
completed if we use Eq.® in the section of free propa- 
gation, and for the action of thin lens wc multiply the 
wavefunction at the plane of the lens by the phase factor 
exp(— ifc|spp/2/), where / is the focal length of the lens. 
It is important to maintain the same dimensionless pa- 
rameter s along the entire optical system in order to use 
the additivity property: Ta+ii = Ta ° J'f3 of the FRFT. 
This is the mechanism behind the implementation of a 
FRFT between two planar surfaces with the optical sys- 
tems reported in 3J|, where (frip) = ^a[^eip)]- In the 



experiment reported in section [TVl we perform FRFT's 
using the "type I" symmetric lens system configuration, 
which is illustrated in FIG.[T]b). This FRFT system was 
orig inally reported in [3^ and is also discussed in detail in 
[3(j | . This FRFT system consists of a lens of focal length 
/ placed symmetrically between the input and output 
planes, at a distance Za from each. One can apply ei- 
ther Fourier optics or geometric optics to verify that this 
system corresponds to a FRFT. Specifically, it is neces- 
sary to define the fractional focal length /' = / sin a, and 
impose that the focal length / and the distance of prop- 
agation Za before and after the lens are related to the 
order a of the FRFT via the relation Za = 2/sin^(a/2). 
The dimensionless position and mom entum coordinates 
for this kind of system are p = \/k/ f'p and q — ^ f /k(\, 
where p, q are the dimensional variables. In operator for- 
malism, the evolution with a FRFT is associated with the 
fractional Fourier operator defined as [soj 



exp —la 



(6) 
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FIG. 1: a) In free diffraction, the field on the curved surface 
with radius Rr can be described as the FRFT of the field on 
surface of radius Re with properly scaled coordinates, b) The 
FRFT can be implemented with a simple lens symmetrical 
system [s^j. 



where p and q arc the dimensionless position and mo- 
mentum operators. This operator is equivalent to the 
evolution operator of the quantum harmonic oscillator, 
with the hamiltonian defined as H = (p^ + q^)/2. Under 
the FRFT operator of order a, the single photon state 
I (/Jo) evolves to \ipa) = Ta\'^o)- The FRFT of the wave 
function ^aip) = {p\^a) is then given by [s^l 



.{P)= I dp' 



P\^a \p') fo{p'), 



(7) 



where the kernel is 



(pI-^oIp') 



.cot a ,2 



X exp I 



cot a 



exp 



-^), (8) 



vector operators evolved under the action of FRFT are 



Pa 



COS a sm ol 
— sin OL cos OL 



(9) 

which illustrates the fact that J- a corresponds to rotation 
of angle a. in phase space [sol . [3^ . 



III. THE ENTANGLED TWO PHOTON STATE 

Let us consider now a pure two photon state whose 
wavefunction in dimensionless coordinates is of the form 



^(Pi, P2) = /(Pi + P2)g{pi - P2) 



(10) 



This state is generally correlated, provided that /(p) and 
g{p) are not identical Gaussian functions. Here it is as- 
sumed that / and g are normalized with respect to pi, p2- 
The state pUj) can be readily produced in a number of 
physical processes [ll|, . It is representative of the two 
photon state at the face of the SPDC crystal, for exam- 
ple, provided that the pump and down-converted fields 
are polarized and nearly monochromatic [4l| . In this case 
/ is given by the spatial profile of the pump field and g 
is the Fourier transform of the phase matching function 
G'(q) = y/2L/KT:H\nc{L\q^/AK) [i^, where K is the 
wave number of the pump beam. In many experimental 
situations, G(q) and g{p) can be approximated by Gaus- 
sian functions. In this case, assuming that the pump laser 
has a Gaussian profile, the position space wave function 
takes the form 



'J'(P1,P2) = 



1 



|P1 +P2 I 



(11) 



Eq. (fTTjl describes the field at the crystal face. In trans- 
verse wave- vector space the wave function is 



*(qi,q2) 



e 



-|qi 



-|qi-q2 



(12) 



which is obtained by taking the Fourier transform of the 
wave function ()lip . Now let us suppose that (t_ ^ (t+, 
so that the photons exhibit a position correlation and a 
momentum anticorrelation. This is indeed what is gen- 
erally produced in SPDC, in which it is not unusual to 
have cr_ ~ cr^/lOO. 



A. Propagation of transverse correlations 



for < |a| < TT. Here = — iexp(jQ!/2)/(27r| sina|). 
Taking the limit a — > (or a — > 27r), one can show 
that {p\Ta\p') = S(p — p') and similarly {p\!Fa\p') = 
S{p + p') for a ^ ±7r [lol. When a = 7r/2, the FRFT 
reduces to the common Fourier transform. When a does 
not lie in the interval Q < \a\ < n, Eq. ([5]) accurately 
represents the FRFT kernel provided one replaces a with 
its value modulo 27r. The transverse position and wave- 



As discussed above, propagation of the down-converted 
fields can generally be described by a FRFT operation. 
Let us assume that photon 1 propagates according to 
an a-order FRFT along axis zi, and photon 2 according 
to a /3-order FRFT along axis Z2- The state \^) after 
propagation is given by 

|*a,/3)=.FW®.Ff (13) 



5 



The two-photon wave function then becomes 
"^aAPi^P^) = (Pi,P2|*Q,/3>, where 



*a,/3(pi,P2) = jj dp'^dp'^ {pl\J'a \p'i) X 

(P2|^/3|P'2)*(P'1,P2), (14) 

and the kernels are defined in Eq. ([8]). To get a sense of 
the action of the FRFT's, let us consider the limiting case 
of an Einstein-Podolsky-Rosen (EPR) state, for which 
f{p) ~ constant and g{p) ~ 5{p), giving *(pi,p2) = 
5{pi — p2)- This situation is approximated by the state 
produced by SPDC when the pump beam can be treated 
as a plane wave. The EPR state is 

1*^™) = // dp,dp25{p, - P2) |P2>2 , (15) 

which presents a perfect correlation, since detection of 
photon 2 at position p projects photon 1 onto a position 
eigenstate \p). After FRFT's, the wave function is 



*a™(pi,P2) =A,A^exp 



. cot ap\ + cot (3p\ 



, cotap^\ /.cot/3p^ 
apexp i — exp i — 



X exp 



ip- 



Pi 

sin a 



P2 

sin/3 



(16) 



Performing the integral, we have 

.cot (a + P) 



*a™(Pl,P2) =A„+/3exp 



'A + pi) 



exp 



; Pi ■ P2 

sin(a + /3) 



(17) 



which is the kernel of an FRET of order a + (3 corre- 
sponding to propagation from an input plane (e.g. pi) 



to an output plane (e.g. P2). The state 



^'5?™^ is then 



exp 



dpidp2 exp 

. Pl ■ P2 



'. C0t(Q + /3) ( 2 , 2\ 
* 7, [P1+P2) 



sin (a + (3) 



\pi)i m) 



(18) 



Using the definition of the FRET kernel ([8]), we note 
that whenever a-f/J = (mod27r), the original state (fT51) 
is recovered. That is, the EPR state ([15]) is an eigenstate 
of operators of the type TaT2-K-cn TaT^-n-a, etc. When 
a-|-/3 = TT (mod27r), the correlated EPR state evolves 
to an anticorrelated EPR state 



1*^™) = // dpidp2^(pi +P2) |pi)i IP2; 



2 • 



(19) 



In this case the detection of photon 2 at p projects pho- 
ton 1 onto the state |— p). Given any propagation char- 
acterized by an FRET Ta on photon 1, one can find a 



transformation Tfi on photon 2 such that a correlation or 
anticorrelation is recovered. Whena-|-/3 = 7r/2 (mod27r), 
this state becomes 



dp\p)i |q(p))2 



(20) 



where |q(p)) c>c Jexp(iq(p) • p) |p) is the momentum 
eigenstate conjugate to |p). State (|20[) presents no in- 
tensity correlation. An equivalent result is found for 
a + P = 3tt/2 (mod27r). We note that the conditions for 
correlation, anti-correlation, and no-correlation depend 
on the sum of the FRET angles of the down-converted 
fields, and not the individual angles a and p. 

This simple picture drawn for the ideal EPR-state 
is followed approximately by the two-photon state in 
Eq. pTj) . For simplicity, let us use the fact that the two- 
photon wave function is factorable in x and y variables: 
*(pi,P2) = ^{Pxt,Px-2)S.{Pyi,Py2)- Then we can consider 
one spatial dimension p for each down-converted field. 
Figures[2]and[3]show the initial state ([TT]) propagated un- 
der different FRFT's using a+ = 4.076 and cr_ = 0.067. 
FIG.Oshows two examples of strong correlations between 
photon 1 and 2 when the ERFT orders satisfy the condi- 
tion a + /3 = (mod27r) and two examples of strong an- 
ticorrelations when the condition \s a + (3 = tt (mod27r). 
For the condition a + [3 = 37r/2, FIG. [3] shows a signif- 
icant decrease of intensity correlations, although in gen- 
eral they do not completely dissapear as is the case shown 
in Eq. ij^n)) for the ideal EPR state. In fact, analytical 
calculation shows that, in order to have no intensity cor- 
relation, i.e. |(pi,P2|*Q,/3)P = /i(pi)/2(P2), the exact 
relation between a and (3 is: 



tan a tan (3 



^2 2 

<7_<T I . 



(21) 



Eq. ([21]) is satisfied by FRET orders such that a + /? = 
7r/2 (mod27r) or a + (3 = 37r/2 (mod27r) only when 
CT_ = 1/(7+. Nevertheless, the intensity correlations 
present in the state Eq. (jlip propagate in a fashion sim- 
ilar to idealized case of the EPR state. 

In the laboratory, one has access to the joint detec- 
tion probability, which in the case of a two-photon state 
corresponds to the fourth-order correlation function [43| 

-Pa,^(pi,p2) = (*„,/3|a^(pi)at(p2)a(pi)a(p2)|«'a,/3) 
= \^o.APuP2)\\ (22) 

and is proportional to the number of coincidence counts 
Ca^pipi, p2)- The conditional probability can be ob- 
tained by the relation 



PaAP2\Pl 



Pa,l3iPl,p2) 

Ppipl) ' 



(23) 



where Pf3{pi) is proportional to the number of sin- 
gle counts Cp{pi). Thus, the conditional probability 
Pa,p{p2\Pi) is also proportional to the number of two- 
photon coincidence counts Ca,i3{pi, P2)- 
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FIG. 4: (color online) Experimental Setup. The distances of 
the FRFT lens systems are di = 42.63cm and d2 — 7.33cm. 
All lenses have focal length / = 25cm. Moveable slits (not 
shown) are placed in front of each detector. 



FIG. 2: Density plot of the joint detection probability calcu- 
lated for the initial gaussian state in Ea. (|lip evolved with the 
following FRFT orders: a) a = 37r/4, P = 57r/4; b) a = vr, 
/3 = tt; c) a = 7r/4, /? = 37r/4 and d) a = 7r/2, /3 = tt/2. A 
strong correlation is present when a+/3 = (mod27r) [plots a) 
and b)] and a strong anti-correlation when a + /3 — n (mod27r) 
[plots c) and d)]. 
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FIG. 3: Density plot of the joint detection probability calcu- 
lated for the initial gaussian state in Eq. (|lip evolved with the 
following FRFT orders: a) a = 7r/4, /3 = 57r/4; b) a = 7r/2, 
f} — TT and c) q: = 37r/4, /3 = 37r/4. In this case, with 
a + P = 37r/2 (mod27r), we have a strong decrease in the 
correlations. 



IV. EXPERIMENT 

We investigated the propagation of EPR-like correla- 
tions experimentally by implementing several FRFT's on 
pairs of entangled photons and registering the coinci- 
dence counts while scanning one of the detectors. The 
experimental setup is shown in Fig. 3) Degenerate twin 
photons with A — 810nm arc generated by pumping a 
5mm long lithium iodate crystal (LilOs) with a lOmW 
c.w. diode laser centered at Xp ~ 405nm. The transverse 
waist of the beam at the laser output was measured to 
be 0.31 ± 0.01 mm. To increase the spatial correlations, 
the beam width is expanded three times using two con- 
focal lenses. The down-converted photons are detected 
by APD photodetectors equipped with lOnm bandwidth 
interference filters centered at 810nm. Moveable hori- 



zontal slits (100/imx3mm) arc placed directly in front of 
each detector in order to scan the vertical position. The 
FRFT's are performed on both down-converted fields us- 
ing the "type I" symmetric lens system configuration re- 
ported in [3J| and shown in FIG. [T]b). The dimension- 
less position and momentum coordinates for this kind 
of system are p = \/k/ f'p and q = \fj^~f^^^ where 
/' ~ 2b/y/2 cm (see below) is the scaled focal length and 
p, q are the dimensional variables. 

Initially, correlation measurements for the near field 
(a = /3 = tt) and far-field (a = /? ~ 7i'/2) correlations 
were obained by fixing one detector at p = and scan- 
ning the other [23|. These correlations correspond 
to the usual position and wave vector variables in the 
source plane. The near field correlations were obtained 
by imaging the exit surface of the crystal on the plane 
of the detectors with 4/ lens systems. For the far-field 
measurements, the usual optical Fourier transform sys- 
tem was used. A sample of the coincidence counts are 
shown in FIG. [51 as a function of the dimensionless vari- 
able p2 ■ The conditional variances are listed in table [D 
Using these results, we can evaluate the EPR inequality 

A2 jpi|p2)A2_Jgi|g2) = 0.20 ±0.01 < \ (24) 
l^l^MPi)^l,Mqi) = 0.14 ± 0.01 < i, (25) 

which shows that the state displays non-local correla- 
tions. Also shown in FIG. [5] are the results using differ- 
ent lens configurations, which give the weakly correlated 
intensity distributions. For example, the tt — 7r/2 distri- 
bution is more than 10 times larger than the tt — tt and 
7r/2 — 7r/2 distributions. 

To evaluate the strength of these correlations under 
different FRFTs, a series of measurements were per- 
formed with various FRFT lens systems. All lenses used 
in the experimental setup have the same focal length 
/ = 25cm. We chose FRFT's with orders a = {^, ^} 
and = {f , where a and (3 correspond to photons 
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FIG. 5: (Color online) Coincidence counts Ca,i3{pi, P2) as a 
function of dimensionless p2 for a.) a = n, [3 = n (red circles) 
and a — 7r/2,/3 = vr (blue squares) and h) a = 7r/2,/3 = n/2 
(red circles) and a — tt, l3 — tt/2 (blue squares). In all cases 
detector 1 is fixed at pi — 0. 



1 and 2. respectively. These FRFT orders sum to either 
TT, 37r/2, or 27r. This choice of angles is especially conve- 
nient, as it maintains /' — 25/-\/2 cm the same for all of 
the FRFT systems used. This is advantageous for several 
reasons: i) to respect the condition of additivity of two 
consecutive FRFT's systems and ii) to use the same scal- 
ing factor for signal and idler fields which is necessary in 
order to describe the FRFT mathematically as a rotation 
in phase space. The scaling parameter for our system is 
y/k/f = 6.62mm-i. 

The various lenses used to implement these FRFTs are 
shown in FIG. ID Three additive FRFT lens systems were 
used to perform the 57r/4 order FRFT. Lens Li is used 
to perform a 37r/4 order FRFT of the field from the exit 
face of the crystal to position 2c?i. Lenses L2 and L3 
each perform a 7r/4 order FRFT, the first from z = 2di 
to z = 2di + 2d2 and the second from z = 2di + 2d2 to 
z = 2di +Ad2- The field at the plane z — 2di +Ad2 is the 
FRFT of order 57r/4 of the field at the exit face of the 
crystal. Lens L5 is used to perform a ^ FRFT, and L4 
was used to perform a 7r/4 order FRFT. By choosing dif- 









'-^ci,f3yP2\h'l ) 


'-^a,l3\Pl\h'2 ) 




2 ) P — 


TT 

2 




1 7 + n 09 




- TV R — 






1 39 + Ofi 








12.3 ± 1.5 




a = 


= 7r,/3 = 


7r 
2 


13.3 ± 2.1 




a = 


4 ' P 


5ir 
4 


0.29 ± 0.01 


0.28 ± 0.02 


a = 


- 3 = 


3ir 

4 


0.21 ± 0.01 


0.31 ± 0.02 


a = 


f,/3 = 


5lT 

4 


9.2 ± 0.7 


13.8 ± 1.2 


a = 




37r 
4 


12.1 ± 1.1 


17.7 ± 1.1 



TABLE I: Conditional variances for all measurement results 
for different FRFT orders a and /3. Variances were obtained 
from gaussian curve fits. 



ferent detector positions and combinations of lenses, we 
could implement several different FRFT's on each down- 
converted field. 

A sample of the experimental results are shown in 
FIG. [HI which displays coincidence counts Ga,f}{pi, P2) 
as a function of the dimensionless coordinate p2- In all 
of the plots the slit of detector 1 is fixed at the ori- 
gin (pi=0). These figures correspond to vertical cross- 
sections along the line pi = of the theoretical density 
plots in FIG. 's [2] and [31 One can see that for the cases 
a + [3 ~ TT, 27r, a narrow coincidence distribution is ob- 
served, indicating either an intensity correlation or anti- 
correlation. When a + f3 = 37r/2, the coincidence pro- 
file is much larger, indicating a much weaker correlation. 
Using Eq. (p3)) . the conditional variances A^^(p2|pi) 
were determined through gaussian curve fits of the coin- 
cidence distributions. Similar measurements and analysis 
were conducted by scanning pi and fixing detector 2 at 
P2 = 0. The dimensionless variances for all results ob- 
tained are presented in table |T1 We note that the variance 
for the weakly correlated distributions are about 10-50 
times larger than the correlated and anti-correlated dis- 
tributions. 

The coincidence distribution 3^{pi\p2) in fact cor- 
responds to the the transverse wave-vector distribution 
C^ !^{qi\q2), since the a = 7r/4 (/? = 37r/4) FRFT dif- 
fers from the a = 3tt/A {(3 = 5tt/A) FRFT by a Fourier 
transform. Thus, with the experimental results shown in 
Fig. [6l we can calculate the EPR inequality ([T]): 



AL 5^(pi|p2)AL 5^(<7i|g2) = 0.0391 < ], (26) 

4 ' 4 4 ' 4 4 



indicating EPR nonlocality. Similiarly, the conditional 
variances A^_^(pi|p2) give 

AL 5^(P2|pi)AL s^felgi) = 0.0352 < ]. (27) 

4 ' 4 4 ' 4 4 

It is clear that the EPR intensity correlation is lost when 
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FIG. 6: (color online) Coincidence counts Ca,p{pi, P2) as a 
function of dimensionless p2 for a) a = 7r/4, /3 — 37r/4 (red 
circles) and a = 37r/4, /3 — 37r/4 (blue squares) and h) a — 
37r/4, /3 = 57r/4 (red circles) and a = 7r/4, /3 = 57r/4 (blue 
squares). In all cases detector 1 is fixed at pi — 0. 



a + (3 = 37r/2, since 

:(Pl|P2)A 



2 ^„ i„ ^A2 

37r 3ir 



,(gi|(Z2) = 244±26> 



A2.^(p2|pi)A^.Jq2|gi) = lll±13> 



(28) 
(29) 



The results summarized in table U show the strength, 
but do not indicate the type of correlation. To investi- 
gate the type of spatial correlation in intermediate FRFT 
planes, we first used lens configurations with FRFT or- 
ders a = ^ and /? = satisfying a + (3 — 2tt. Exper- 
imental results are shown in FIG. [7] a). Initially, the slit 
in front of detector 1 was placed at the origin {pi ~ 0) 
and the slit in front of detector 2 was scanned vertically. 
The measured coincidence counts are plotted in black tri- 
angles in figure [7] a) and the maximum of the gaussian 
fit is at pi = ± 0.03. We then displaced one of the shts 
by 300 ± 5/im, which corresponds to a dimensionless dis- 
placement of pi = 1.99 ± 0.03. Coincidence counts were 
again measured while the slit of detector 2 was scanned. 
Coincidence counts are in blue squares in FIG. [7] a) and 
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FIG. 7: Coincidence counts Ca,0[pi, P2) as a function of di- 
mensionless p2 for a) Q = 37r/4 and /3 = 57r/4 and and b) 
a — 7r/4 and /3 — 37r/4. In both figures, the black triangles 
correspond to pi = 0, the red circles to pi = —1.99 ± 0.03 
and the blue squares to pi = 1.99 ± 0.03. Figure a) thus 
shows a correlation, while b) shows an anti-correlation of the 
detection positions. 



the maximum of the gaussian fit is at p2 ~ 1.93 ± 0.03. 
Slit 1 was then moved -300±5^m, (pi = -1.99±0.03), 
and slit 2 was scanned. The maximum of the coinci- 
dence counts occurred at p2 = —1.94 ± 0.03. We thus 
observe a strong correlation between the transverse coor- 
dinates for this configuration satisfying a + (3 = 2tt. The 
same procedure was performed for the lens configuration 
{a = = f}, which satisfies the anti-correlation 

condition a + (3 — n. The results are shown in FIG. [7| 
b). We observe similar displacement of the coincidence 
peaks, however in this case the maxima of the gaussian 
fits are anti-correlated with the position of the slit of the 
other detector. 



V. CONCLUSION 

We have used the fractional Fourier transform to study 
the propagation of the transverse intensity correlations of 
the two-photon state produced from parametric down- 
conversion. The transforms were implemented with sim- 
ple lens systems. Our theoretical and experimental re- 
sults show that the propagation of the transverse cor- 
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relations of highly-correlated two-photon states depends 
upon the the sum of the transform orders of the down 
converted fields. For a + (3 = (mod27r), the origi- 
nal intensity correlation at the source is recovered, while 
for a + P ~ TT (mod27r), an intensity anti-correlation 
is observed. For a + /3 = tt/2 (mod27r) or a -I- /? = 
3tt/2 (mod27r), almost no correlation is present. Ana- 
lytical results were obtained for the propagation of the 
ideal EPR state, and numerical calculations along with 
our experimental results show that the down-converted 
photons display a similar behavior. The EPR correlation 
present in the two-photon state was confirmed for sev- 
eral different orders of the fractional Fourier transforms 
through violation of an inequality. These results apply 



to spatially correlated photons obtained from any source, 
as well as correlations present in other physical systems, 
and should be useful for engineering spatial correlations, 
as well as fundamental studies of quantum nonlocality 
and entanglement. 
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